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Abstract 

Following the newly formulated notion of form invariance of the neutrino mass 
matrix, a complete model of leptons is constructed. It is based on a specific unitary 
3x3 matrix U in family space, such that U 2 is the simple discrete symmetry v e — > 
— z^ e , i^j <-> i/ T . Thus £/ also generates the cyclic group Z±. The charged-lepton mass 
matrix is nearly diagonal while the neutrino mass matrix is of the form suitable for 
explaining maximal (large) mixing in atmospheric (solar) neutrino oscillations in the 
context of three nearly degenerate neutrino masses. Observable lepton flavor violation 
is predicted. Quarks may be treated in the same way as the charged leptons. 



To understand the form of the neutrino mass matrix a new idea has recently been 
proposed [T|. It is postulated that there is a specific 3x3 unitary matrix U with U n = 1 
such that 

UM U U T = M u , (1) 

and for some N < n, the matrix U represents a well-defined discrete symmetry in the v e ,n,r 
basis, for which the charged-lepton mass matrix A4i is diagonal. However, since each neu- 
trino belongs to an SU(2)l x U(1)y doublet, the corresponding left-handed charged leptons 
(e, /i, r) must also transform under U. A complete theory must then reconcile the apparent 
contradictory requirement that M. v satisfies Eq. (1), but Aii does not. The resolution of 
this conundrum is in the soft and spontaneous breaking of the symmetry supported by U, 
as was done for example in the A 4 model [3 El of degenerate neutrino masses. In this paper 

we show how this may be achieved with the specific unitary matrix 

/ i/V2 -i/V2\ 



U 



i/y/2 1/2 1/2 
\-i/V2 1/2 1/2 J 



(2) 



where 



U 



f-l 0\ 
1 

v o i oy 



(3) 



i.e. the simple discrete symmetry 



-Z/ e , Uf, U T 



(4) 



The matrix U of Eq. (2) is thus a "square root" of this discrete symmetry. We see immediately 
also that U 4 = 1; hence our model is a specific realization of the cyclic group Z 4 as a family 
symmetry. 



The most general Majorana mass matrix A4 U is of the form 

(A D E\ 
M u = D B F . 
K E F Cj 



(5) 



Using Eqs. (1) and (2), we find that it becomes 

^ 

(6) 



M, 



(AO \ 

B A+B . 
v A + B B j 

This mass matrix has eigenvalues A, —A, and A + 2B, corresponding to the eigenstates u e , 
(v„ — Vr)/y2, and + u T )/y/2. We see immediately that — v T mixing is maximal with 



I'- 



Am 2 atm = (A + 2Bf - A 2 = AB{A + B), (7) 
which is suitable for explaining atmospheric neutrino oscillations [3]. If A << B, we have 



the hierachical structure of neutrino masses and B = y Am 2 . tm /2 ~ 0.025 eV. On the other 
hand, if B << A, we have the more interesting scenario of three nearly degenerate neutrino 
masses, with the prediction that A is large enough to be measured by neutrinoless double beta 
decay. As for solar neutrino oscillations [3], we can obtain the large- mixing- angle solution 
by invoking flavor-changing radiative corrections Details will be presented in a later 

paragraph. 

The leptonic Yukawa couplings of our model is given by 

+ f^-v^lt + H.c, (8) 

where we have adopted the convention that all fermion fields are left-handed, with their right- 
handed counterparts denoted by the corresponding (left-handed) charge-conjugate fields. We 
have also extended the Higgs sector of the Standard Model of particle interactions to include 
three doublets ((f)®, (f)j) and one very heavy triplet The electroweak SU(2)i x 

U(l)y gauge symmetry will be spontaneously broken mainly by the vacuum expectation 
value of one particular Higgs doublet. As a result, the vacuum expectation values of the 
other two Higgs doublets can be naturally small jHj and that of the triplet even smaller [Zj. 



We now assume that Ly is invariant under the transformation 



I 



U, 



>3i 



hi 

(f 



This means 



U T hU = h, U T fU = f 



T fk 



(9) 
(10) 

(11) 



resulting in 

(a \ / a k d k -d k \ 

h= b a + b , f k = -d k b k a k + b k 
V a + b b J \ d k a k + b k b k ) 

Note that h has no d terms because it has to be symmetric. 



(12) 



The neutrino mass matrix Ai u is then given by 

M u =2h(e), 



(13) 



whereas the charged-lepton mass matrix A4i linking li to 1% is 

/ a x v x + di(v 2 - v 3 ) a 2 v 1 + d 2 {v 2 - v 3 ) a 3 Vi + d 3 (v 2 - v 3 ) \ 

-d\Vx + h{v 2 + v 3 ) + aiv 3 -d 2 v x + b 2 (v 2 + v 3 ) + a 2 v 3 -d 3 Vi + b 3 (v 2 + v 3 ) + a 3 v 3 , 

\ d x vx + a x v 2 + b x (v 2 + v 3 ) d 2 Vi + a 2 v 2 + b 2 (v 2 + v 3 ) d 3 v t + a 3 v 2 + b 3 (v 2 + v 3 ) J 

(14) 

where Vi = (0°). Assume dk, b k « a k and Vi i3 « v 2 , then all elements in the first, second, 
and third rows are of order dv 2 + avi, bv 2 + av 3 , and av 2 respectively. It is clear that they may 
be chosen to be of order m e , m M , and m T , in which case Aii will become nearly diagonal by 
simply redefining the l c k basis. The mixing matrix Vl in the li basis (such that VlM.iA4.IvI 
is diagonal) will be very close to the identity matrix with off-diagonal terms of order m e /m M , 
m e /m T , and m M /m r . This construction allows us then to consider Ai u to be in the basis 
(z/ e , z/ M , z/ r ) to a very good approximation. 

We now understand why it is sensible pQ| to consider Eq. (1) as a condition on Ai u . The 
key lies in the fact that M. v comes from neutrino couplings to a single field £° which is 



which shows clearly that neutrino masses may be of order 1 eV or less if M?/fi ~ 10 13 



invariant under U, whereas Aii comes from couplings to 0? 2 3 which are not. Let the scalar 
trilinear coupling of £ to $2 $2 be /i, then [Zj 

1 

GeV. Similarly, 1^3 can be small compared to t>2 if M 2 < but the Mf 3 terms in the Higgs 
potential are large and positive, in which case [U] 

^^T' ^Hup (16) 

where [i\ 2 and //| 3 are the coefficients of the ${$2 and $2^3 terms respectively. 

Going back to A4 U of Eq. (6), we now consider how solar neutrino oscillations may arise 
in the 2x2 submatrix spanning v e and {y^ — v T ) / \/2, i.e. 

Mo -J- 

Consider the most general radiative corrections to A4, i.e. 

R=h r "), (18) 
VI2 r 22 J 



then M. becomes 



( i + s)M(i + r ) ^( 1 , +2r " I 2 TV < i9 > 

VrJ 2 — r 12 -l-2r 22 / 
In terms of the full radiative correction matrix 

rn = r ee , r 22 = ^(r M(U + r rr ) - Re^^), (20) 

r* 2 - r i2 = -2ilm{ri 2 ) = -iV2[Im(r efl ) - Im{r eT )\. (21) 

Thus the radiatively corrected M. has eigenvalues 

m 1>2 =A(l + r n + r 22 =F y/ (r n - r 22 ) 2 + 4[/m(r 12 )] 2 ) (22) 



corresponding to the eigenvectors v e cos 8—i(u^ — u T ) sin9/\/2 and v e cos 6+i{u^ — v T ) cos 6/ \[2 
respectively, where 



„ r n -r 22 + J(r n - r 22 ) 2 + 4[/m(r 12 )] 2 

tan# = v , T . — — — , 23 

2|/m(r 12 )| ' 1 ; 

with r 22 — r n > 0. Since tan 2 9 ~ 0.46 is desirable jH] for understanding solar neutrino 
oscillations [5 , flavor changing (r 12 ^ 0) and flavor nonuniversal (r n 7^ r 22 ) interactions are 
required. Specific examples have already been proposed [HIS]- As for Am 2 oZ , it is given here 
by 



Ami, = m 2 2 - m\ ~ 4A 2 ^ (r n - r 22 ) 2 + 4[/m(r 12 )] 2 . (24) 

For |r| ~ 10~ 3 which is a typical size for radiative corrections and Am^ ~ 6.9 x eV 2 
jH], \A\ ~ 0.13 eV is then obtained. [The recent WMAP result implies 0| an upper bound 
of 0.23 eV on \A\ from cosmological considerations.] Given that \A\ is also the effective 
neutrino mass measured in neutrinoless double beta decay with a present upper bound of 
about 0.4 eV, this is an encouraging prediction for future experiments ^U]- If A ~ 0.13 eV, 
then using Eq. (7), we find B ~ 0.0048 eV and B/A ~ 0.037 which is of the same order as 
fn^/rrir ~ 0.059, as suggested by M.\ of Eq. (14). 

With flavor changing radiative corrections, the U e s entry of the neutrino mixing matrix 
becomes nonzero. It is given here by 

Ue3 ~ ' (25) 

which may be as large as the experimental upper bound [TT] of 0.16. However, it is real so 
that there is no CP violation. This is in contrast to the case of the A^ model [B], where C/ e 3 
has to be purely imaginary [T2] . 

Going back to the Yukawa couplings of the leptons to the 3 Higgs doublets given by 
Eq. (12) and assuming the hierarchy dk << bk << at and taking the limit that only t> 2 is 
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nonzero, we have Aii of Eq. (14) simply given by 

/ di d 2 d 3 \ 
Mi ~ v 2 bi b 2 h 
\ai a 2 a 3 / 

whereas $1 and $3 couple to 1^ according to 



(26) 
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a 1 


a 2 


«3 \ 
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—di 


-c/ 2 


-d 3 \ 




—d\ 


-d 2 


-d 3 , 




a\ 




a 3 
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di 


d 2 


d 3 ) 
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61 


&2 


&3 / 



(27) 



respectively. After rotating .M; of Eq. (26) in the lj basis to define the state corresponding 
to r, we see immediately from Eq. (27) that the dominant coupling of $1 is (m T /v 2 )er c and 
that of $3 is (m T / v 2 ) fir c . Other couplings are at most of order m fJi /v 2 in this model, and 
some are only of order m e /v 2 . We thus have a natural understanding of the smallness of 
flavor changing decays in the leptonic sector, even though they should be present and are 
potentially observable. 

Using Eq. (27), we see that the decays r~ — > e~e + e~ and r~ — > e~e + fi~ may pro- 
ceed through cf)\ exchange with coupling strengths of order m^m T /v 2 ~ (g 2 /2)(m tl m T /M- 2 v ), 
whereas the decays t~ — > fi~fi + fi" and r~ — > fi~fi + e~ may proceed through 0° exchange 
also with coupling strengths of the same order. We estimate the order of magnitude of these 
branching fractions to be 



B 



1 B(t — > fiisv) ~ 6.1 x 10 1 



(28) 



V Ml, J v ' V M i,3 / 

which is well below the present experimental upper bound of the order 10~ 6 for all such rare 
decays [T5] . 

The decay /i~ — > e~e + e~ may also proceed through 0? with a coupling strength of order 



rri' l Jv\. Thus 



m, 



eee) ~ — £ ~ 1.2 x 10~ 12 



100 GeV \ 4 
M 1 ) 



(29) 



which is at the level of the present experimental upper bound of 1.0 x 10~ 12 . The decay 
u — > ey may also proceed through 0° exchange (provided that Recj)® and Irrufy have different 
masses) with a coupling of order m^m T /vl. Its branching fraction is given by [2] 



Using the experimental upper bound jH] of 1.2 x 10 , we find M e ff > 164 GeV. 

In the quark sector, if we use the same 3 Higgs doublets for the corresponding Yukawa 
couplings, the resulting up and down mass matrices will be of the same form as Eq. (14). 
Because the quark masses are hierarchical in each sector, we will also have nearly diagonal 
mixing matrices as in the case of the charged leptons. This provides a qualitative under- 
standing in our model of why the charged-current mixing matrix linking up quarks to down 
quarks has small off-diagonal entries. 

Once (f)\ or <f)® is produced, its dominant decay will be to r ± e =F or r ± yU =F if each couples 
only to leptons. If they also couple to quarks (and are sufficiently heavy), then the dominant 
decay products will be tu or tc together with their conjugates. As for 02, h will behave very 
much as the single Higgs doublet of the Standard Model, with mostly diagonal couplings to 
fermions. 

In conclusion we have constructed a complete theory of leptons where the neutrino mass 
matrix M. v may be derived from the requirement that UM. V U T = Ai u , where U is a specific 
3x3 unitary matrix in family space such that U 2 is the simple discrete symmetry u e — > —u e , 
v T . We obtain three nearly degenerate neutrino masses with maximal mixing for 
atmospheric neutrino oscillations. Solar neutrino oscillations are induced by flavor-changing 
radiative corrections. As a result, neutrino less double beta decay is predicted to occur at the 



B(u -> ej) 



3a m\ 



(30) 



where 




(31) 
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0.1 eV range. There are also three Higgs doublets in this model, two of which have dominant 
flavor- changing couplings proportional to m T and may be easily observed at future colliders. 
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